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Abstract 

The problem of maximum rate achievable with analog network coding for a unicast communication over a layered relay 

network with directed links is considered. A relay node performing analog network coding scales and forwards the signals 

received at its input. Recently this problem has been considered under certain assumptions on per node scaling factor and received 

Cn SNR. Previously, we established a result that allows us to characterize the optimal performance of analog network coding in 

^~~^ network scenarios beyond those that can be analyzed using the approaches based on such assumptions. 

^^ The key contribution of this work is a scheme to greedily compute a lower bound to the optimal rate achievable with analog 

Cn network coding in the general layered networks. This scheme allows for exact computation of the optimal achievable rates in a 

5_j wider class of layered networks than those that can be addressed using existing approaches. For the specific case of Gaussian 

^ I A^-relay diamond network, to the best of our knowledge, the proposed scheme provides the first exact characterization of the 

.^f^ optimal rate achievable with analog network coding. Further, for general layered networks, our scheme allows us to compute 

^ optimal rates within a constant gap from the cut-set upper bound asymptotically in the source power. 

m 

psj I. Introduction 

^^ Analog network coding (ANC) extends to multihop wireless networks the idea of linear network coding 1 1 1, where an inter- 

^_( mediate node sends out a linear combination of its incoming packets. In a wireless network, signals transmitted simultaneously 
I— ( by multiple sources add in the air. Each node receives a noisy sum of these signals, i.e. a linear combination of the received 
C/3 signals and noise. A communication scheme wherein each relay node merely amplifies and forwards this noisy sum is referred 
I ^1 to as analog network coding [Tl, f3]. 

The rates achievable with ANC in layered relay networks is analyzed in p), B). In p), the achievable rate is computed 
C^ under two assumptions: (A) each relay node scales the received signal to the maximum extent subject to its transmit power 
>^ constraint, (B) the nodes in all L layers operate in the high-SNR regime, where uiinkei PR,k > l/i5, 1 = \, . . . ,L for (5 > 0, 
where PR^k is the received signal power at the /c* node. It is shown that the rate achieved under these two assumptions 
psi approaches network capacity as the source power increases. The authors in H) extend this result to the scenarios where the 
^Sl nodes in at most one layer do not satisfy these assumptions and show that achievable rates in such scenarios still approach the 
-^ network capacity as the source power increasesM 

(^^ However, requiring each relay node to amplify its received signal to the upper bound of its transmit power constraint results 

C^ in suboptimal end-to-end performance of analog network coding, as we show in |6|, |7|. Further, even in low-SNR regimes 
^""^ amplify-and-forward relaying can be capacity-achieving relay strategy in some scenarios, |^5|. 

^ In this paper we are concerned with analyzing the performance of analog network coding in general layered networks, 

k> without the above two assumptions on input signal scaling factors and received SNRs. However, such a characterization of the 

rS performance of analog network coding results in a computationally intractable problem in general [41, 16 J. 

C^ In I?), we establish that a globally optimal set of scaling factors for each node, i.e. a choice of relaying strategies that 

optimizes end-to-end throughput over all ANC strategies, can be computed in a layer-by-layer manner. This result allows us 

to computationally efficiently characterize exactly the optimal ANC rate in a large class of layered networks that cannot be 

addressed using existing approaches under the assumptions A and B. Further, for general layered relay networks, this result 

significantly reduces the computational complexity of computing a set of non-trivial achievable rates. 

However, even layer-by-layer computation of a network-wide scaling vector that maximizes the end-to-end ANC rate for 
general layered networks is a computationally hard problem. In this paper, we propose a greedy scheme to bound from below 
the optimal rate achievable with analog network coding in general layered networks. The proposed scheme allows us to exactly 
compute the optimal ANC rate in a much wider class of layered networks than those that can be so addressed using existing 
approaches, including our approach in |[7l. In particular, for the Gaussian iV-relay diamond network fSl, the proposed scheme 
allows us to exactly compute the optimal rate achievable with analog network coding. To the best of our knowledge, this is 
the first characterization of the optimal ANC rate for Gaussian diamond network. Further, for general layered networks, our 

' However, it is assumed tliat tlie noises at the nodes in this particular layer are independent, resulting in the computed ANC rate overestimating the optimal 
ANC rate in general. 
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Fig. 1. Layered network with 3 relay layers between the source 's' and destination 't'. Each layer contains two relay nodes. 



scheme allows for the computation of the optimal rates within a constant gap from the cut-set upper bound asymptotically in 
the source power. 

Organization: In Section III] we introduce a general wireless layered relay network model and formulate the problem of 
maximum rate achievable with ANC in such a network. Section III addresses the problem of maximum ANC rate achievable 



in a Gaussian A^ -relay diamond network and shows that a greedy scheme optimally solves this problem. In Section IV we first 
generalize the greedy scheme for Gaussian diamond networks to characterize the optimal performance of a specific subnetwork 
of the general layered network. We then construct a scheme to bound from below the optimal performance of ANC in general 
layered networks. Section IV] illustrates that the proposed scheme leads to the exact computation of the maximum ANC rate in a 
specific class of layered networks and tight characterization of the optimal rate in the general layered networks asymptotically 



in the source power. Section VI concludes the paper. 



II. System Model 

Consider a {L + 2)-layer wireless network with directed links. The source s is at the layer '0', the destination t is at the 
layer 'L + 1', and the relay nodes from the set R are arranged in L layers between them. The Z* layer contains ni relay nodes, 






jI=i "-i — ^^- ^^ instance of such a network is given in Figure 111 Each node is assumed to have a single antenna and operate 
in full-duplex mode. 

At instant n, the channel output at node i,i E RU {t}, is 






hjiXj[n\ + Zi[n] 



-oo < n < oo, 



(1) 



where Xj [n] is the channel input of the node j in the neighbor set JV{i) of node i. In ([T]i, hji is a real number representing 
the channel gain along the link from node j to node i. It is assumed to be fixed (for example, as in a single realization of 
a fading process) and known throughout the network. The source symbols Xs[n],—oo < n < oo, are i.i.d. Gaussian random 
variables with zero mean and variance P, that satisfy an average source power constraint, Xs[n] ^ Af{0, Pg)- Further, {zJti]} 
is a sequence (in n) of i.i.d. Gaussian random variables with Zi[n] ^ Af{0, a"^). We also assume that Zi are independent of the 
input signal and of each other We assume that the i* relay's transmit power is constrained as: 

£'[x-[n]] < Pj, -cx) < n < oo (2) 

In analog network coding each relay node amplifies and forwards the noisy signal sum received at its input. More precisely, 
a relay node i at instant n + 1 transmits the scaled version of yi\n], its input at time instant n, as follows 



a;,[n + l]=Ay,W, < ^f < ^^ 



2 

i.max 



P^/P^ 



R,ii 



(3) 



where Pui is the received power at the node i and choice of the scaling factor Pi satisfies the power constraint (|2|. 

One important characteristic of layered networks with unidirectional links is that all paths from the source to destination 
have same number of hops. Also, each path from the z*,i G R, relay node to the destination has the same length. Therefore, 
in a layered network with L layers, all copies of a source signal traveling along different paths arrive at the destination with 
time delay L and all copies of a noise symbol introduced at a node in V''^ layer arrive at the destination with time delay 
L — i + 1. Therefore, the outputs of the source-destination channel are free of intersymbol interference. This simplifies the 
relation between input and output of the source-destination channel and allows us to omit the time-index while denoting the 
input and output signals. 



Using ([T]i and ([3]), the input-output channel between the source and destination can be written as 



yt 



L 



{ii,...,iL)eKL 



1=1 j=l 



J2J2 E /3y/io...../3. 



Li, il-Li, 



{ii,...,iL)£Kij^L 



Zlj + Zt, 



where Kl, is the set of L-tuples of node indices corresponding to all paths from the source to the destination with path delay 
L. Similarly, Kij^^-i+i, is the set of i — / + l-tuples of node indices corresponding to all paths from the j* relay of l'^ layer 
to the destination with path delay L — I + 1. 

We introduce modified channel gains as follows. For all the paths between the source s and the destination t: 

hs ^ E hsi^Pi^hi^i^ . . . Pi^hi^t (4) 

(ii,...,iL)eKL 

For all the paths between the j* relay of /* layer to the destination t with path delay L — / + 1: 

hij — 2_^ PijhijAi ■ ■ ■ hi^hLi^^t (5) 

In terms of these modified channel gains^ the source-destination channel in ^ can be written as: 

L ni 
1=1 3 = 1 

Problem Formulation: For a given network-wide scaling vector (3 = {liu)i<i<L.i<i<ni, the achievable rate for the channel 
in (|6]l with i.i.d. Gaussian input is (BJ, L4J, L6J): 

I{Ps,(3) = {l/2)\og{l + SNRt), (7) 

where SNRt, the signal-to-noise ratio at the destination t is: 

P K^ 
SNRt^^ ^ (8) 

-' 1 + Eti e;li hi 

The maximum information-rate Ianc{Ps) achievable in a given layered network with i.i.d. Gaussian input is defined as the 
maximum of I{Ps.,(3) over all feasible /3, subject to per relay transmit power constraint ([3]). In other words: 

lANciPs)''= max /(P.,/3) (9) 

It should be noted that Pu.max (the maximum value of the scaling factor for i* node in the ^* layer) depends on the scaling 
factors for the nodes in the previous / — 1 layers. 

Given the monotonicity of the log(-) function, we have 

l^iopt — argmax I{Ps,f3) = argmax SNRt (10) 

Therefore in the rest of the paper, we concern ourselves mostly with maximizing the received SNRs. 

In |J7), we discussed the computational complexity of exactly solving the problem (|9| or equivalently the problem ( [TO] i. 
Further, we also introduced a key result |7, Lemma 2] that reduces the computational complexity of the problem of computing 
f3opt by computing it layer-by-layer as a solution of a cascade of subproblems. This result allowed us to characterize the optimal 
end-to-end rate achievable with analog network coding in communication scenarios that cannot be so addressed using previous 
approaches. However, each of these subproblems itself is computationally hard for general network scenarios as it involves 
maximizing the ratio of posynomials |10|, |11|, which is known to be computationally intractable in general |11|. Therefore, 
in this paper, we introduce a greedy scheme that optimally solves these subproblems and consequently the problem ( fTO] ) for 
a large class of layered networks that cannot be addressed with current schemes. For general layered networks, the proposed 
scheme allows us to tightly bound from below the optimal ANC performance. However before discussing this scheme, we 
motivate it by computing the maximum rate of information transfer achievable with analog network coding over the diamond 
network with N relay nodes. 

^Modified channel gains for even a possibly exponential number of paths as in BJ and (3) can be efficiently computed using line-graphs |9|. Further, the 
number of such modified channel gains scales polynomially in the size of the graph being considered. 




Fig. 2. A diamond network with N relay nodes. 



III. Diamond Network: The optimal rate achievable with analog network coding 

Consider the diamond network of Figure l2] We can consider diamond network as a layered network with only one layer of 
relay nodes. Then using Q, (|5]l, and ([S), we compute the SNR at the destination t for any scaling vector /3 = (/3i, . . . , (3^) as 



SNRt 



^' 1 + Eti 



W al-ul 



(11) 



Using (|9|, the problem of computing the maximum ANC rate for this network thus can be formulated as 



max SNRt, (12) 

where Anax = (;9i,max ■ ■ • , /3jv,max) with /Sf^^^ = Pt/{hl,Ps + cr^),i eJV,M={l,...,N}. 

Equating the first-order partial derivatives of the objective function with respect to /3i,i £ Af, to zero, we get the following 
iV + 1 conditions for local extrema: 



y^ hsiPihu = 

/J. = hsi|h^t L^ y^ ^2^2 



(13) 
(14) 



'PiPi ~ g^gfl^' denote the second-order partial derivatives of SNRt with respect to j3i and /3j, i,j £ Af and H{/3) 



Let SNR^ 
denote the determinant of N x N Hessian matrix 

First, consider the set of stationary points Sp = {j3 : (3 satisfies ([T3|)}. For all points in Sp we can prove that 

SNRf,^p^ > 
H{f3) = 

Therefore, the second partial derivative test to determine if the points in Sp are local minimum, maximum, or saddle points 
fails. However, we can establish that for every (3 E Sf3, the following set of conditions holds 

dSNRt 



dSNRt 



aA 



< 0, if ^ hsihitSi < 0, 

> 0, if ^ hsih^tSi > 0, 
Hif3) > 0, if ^ h,,hu6^ < 0, 

Hif3) > 0, if ^ h,,huS^ > 0, 

leAf 



(15) 

(16) 

(17) 
(18) 



for all S — {5i, . . . , S^) -^ 0. In other words, ^T5\ and ([T6]l imply that the slope of the function changes sign at J2ieJ\f ^sihuSi — 
0, and ^n\ and ( fTSJ ) imply that the convexity of the function, however, does not change at J2ieJ\f ^^sihuSi — 0. Therefore, 
together these imply that ( [T3] l leads to a local minimum of the objective function. 



Next, consider the set of points defined by ( [T4] i. For all such points we can prove that 

SNRp^p^ < 
iJ(/3) > 

Therefore, from the second partial derivative test the objective function attains its local maximum at the set of points 
characterized by ( [T4| i above. However, no real solution of the simultaneous system of equations in ( [T4| l exists. In other 
words, no solution of ( [T2| i exists where all relay nodes transmit strictly below their respective transmit power constraints. This 
is illustrated by the following example. 

Example 1 (Three node Diamond Network): Consider the Gaussian diamond network of Figure l2] with three relay nodes. 
For this network, ( [T4] l results in: 

ft' t Hh"'^l' Rh (1 + ft''-!. + feX) (21) 

hsiPihit + hs2P2h2t 

Substituting ( |20l ) in ( |2T] i, after a little algebraic manipulation we get: 

Solving this gives three solutions for fS^, namely: 

/?3,i = ^J^(l + /??/.?,) (22) 

hsiPihit 

n -hs^hsiPihu + lhs2Vih^,2 + fes3)(l + I^Mt) + {hslPlhu^ . ^- 



/i3t(/»L + /^^3) 

Substituting each of (|22]|, (|23|l, and (p4]i in (|20jl results in three corresponding solutions for /32, namely: 



P3,2 - 1, ,1.2 , 1,2 \ ,2-V-i (/4) 



^ _ hnilh2t , ^2.2 N, 

fl-slPl/llt 

Therefore, we have three possible solutions for optimal {^2,^3), i-e. (/32,i, /^s.i), (/?2,2,/33,2), and (/?2,3,/33,3)- 
Substituting (/32,i,/33 1) in ( [T9| ) results in 

o2 "-s2 + "-^3 

which leads to complex valued solutions for optimal /3i. 

Similarly, substituting (/32,2, /?3,2), and (/32,3,/33.3) in ([19]) results in 

02 _ 1 + /»?2 + /»'3 

Pi - i 



hl{hl2 + hl,) 

which also leads to complex valued solutions for optimal /3i. 

This allows us to conclude that no real solution of the system of simultaneous equations in ([T9l)-(|2T|l exists. ■ 

The above discussion implies that all points satisfying ( [T3| l lead to the global minimum of the objective function in ( [T2| ) 
and the global maximum of the objective function occurs at one of the N hyperplanes (of dimension iV — 1) defined by 
Pk = Pk,max, k G Af ■ Ncxt wc identify this hyperplane and characterize the corresponding optimal solution. 

Consider the system of simultaneous equations in ( [T4] l on the {N — 1) -dimensional hyperplane defined by /3k — I3k,max- 



hit iT,skPk,maxi^kt + /^j^J\f\{i,k} '^sjPjhjt 



(25) 



Note that the solution of the above system of equations is the set of scaling-factors for the nodes in the set of relay nodes 
Af\ {k} that maximizes SNRt on hyperplane /3k = Pk.max- Solving the system of equations in ( p5| l results in the following 
set of optimal solutions for j3i on hyperplane (3k = /3fc,; 



■.,max' 



ak "si '^ ~^ h'k,max"'kt ■ . r i r , -, /T<;\ 

^^=^^r—o 7— ,*eA/\{fc} (26) 

*^it '^skPk^max'^kt 

However, the optimal scaling factors in p6l ) for A^ — 1 nodes are computed without considering the upper bound Pi^max on 
each /3i,i € J\f\ {k}. Therefore, taking into consideration the upper bound on the scaling factor for each node, the modified 
solution is computed as per the following lemma. 

Lemma 1: The optimal scaling vector /3q ^ — (/3f ^, . . . , /3^ ^) for N nodes on f3k = (3k,max hyperplane such that each 
scaling factor satisfies the corresponding upper bound on its maximum value is given as 



ok 
^i,opt 



k 



Pi^max T ^ ^ >J 

hsi ^ + 1^'ieS'' Pi,max'^i 



i^S'' 



■ijeS'' r'j,max"'jt 
i^it Z_^j(zgk il'sj Pj.max'^jt 

where S'^ is the set of nodes such that on hyperplane (3^ — /3k,max^ the optimal value of the scaling factor of a node is saturated 
to its corresponding upper bound, S'' — {k} U {i : /3f > (3i^max,i G A/'\ {fc}}- 

Proof: Following the argument similar to the one used to prove the global extrema properties of ( pj) and ( [l4j i, we can 
prove that on the /3k — (3k.max hyperplane, the SNRt achieves its global minimum at a hyperplane defined by 



E 

i£J^\{k} 



/ ^ '^siPi'^i'i 



and its global maximum at the points defined by (3^ given in (|26]l. 

Let M'^ denotes the set of nodes for which (3^ computed in ( |26| l is greater than or equal to the corresponding upper bound 

Pi,max on the maximum value of the scaling factor, i.e. M^ = {i : (3^ > Pi^max: i E J^\ {k}}. For all such f3'^', i e M'^, after 
proving that li''^^ I „ > 0, we set (3^ = f3i max and update S'', the set of nodes such that on hyperplane (3k = (3k max, 

Opi ' Pi, max ' ' 

the optimal value of each node is saturated to its corresponding upper bound; as follows: S^ = S*' VJ M^. As (3^ computed in 
( |26| l for a node i ^ S^ may no longer be optimal after the above re-assignment of /3f , i £ M'^, we need to solve the following 
simultaneous system of iV — \S'^\ = N — |Af'^| — 1 equations with i £ Af\S'': 

, 1 + E Plmaxh% + E ^]h% 

l^it 2-^ f^sj Pj.raax'^jt 'T' /_^ sj Pj '^jt 



Solving this system of equations results in 



jss'' j<^s''u{i} 



/...l + E.e-/^?-"^'^' 



ak _ '^si ~r Z^j^s'' l^j,max jt , „k /^on 

Pi,o-pt — T~"y^ Z o I '* y^ "^ y'^°> 

'''it 2-^i^S^ ''^sj Pj^rnax'^jt 

Some of the recomputed scaling factors (3^ ^,1 ^ S^ may violate the corresponding upper bound on their maximum value. 
All such nodes are added to set S^, thus updating it. Then, the system of equations in ( |27] i is solved again for this updated 
set S^ . This iterative process continues until none of the recomputed (3i in ( |28| l violates its corresponding upper bound. This 
iterative process is presented formally in terms of an algorithm: Algorithm 1, given on the top of the next page. 

Note that Algorithm 1 always halts with either S'' = J\f \ {k} or S'' d M \ {fc} and /?f < (3i^rnax,i €Af\S''. ■ 

Using Lemma fll for each of N hyperplanes, defined as (3^ — (3k.max, k e M, we can compute /3q j, the set of scaling 
factors for all nodes at which SNRt attains its maximum on (3^ — (^k.max hyperplane. Then the hyperplane at which SNRt 
attains its global maximum is identified as follows: 
Proposition 1: The hyperplane at which SNRt attains its global maximum is defined as 

fc* = argmax5iVi?t(/3^ t) 

k£j\f 

Combining Proposition fl] and Lemma fl] we can characterize the scaling vector fiopt that solves the problem ^V2\ as follows. 



Algorithm 1 

1 Initialization: 

5**^ — {k}, the set of nodes whose scaling factors are saturated to their respective upper-bounds 

on hyperplane ^fe = (3k,max- 

U'' ~ J\f\ {k} , the set of nodes whose scaUng factors are not saturated to their respective upper-bounds 

on hyperplane I3k = (3k,max- 

2 Compute /?,f = ^^t%?^l|^!^,zet/^ 

3 while(3/3f >/3,,„a,,zeC/'=) 

4 Compute M'' = {i : pf > f3,,^ax,i e [/'=}. 

6 S"'^ = S"=UM'=. 

7 U'' = U^\ M^. 

8 Compute /3f = ^ 1 + E,-gs>^-^|,^a.^^t ^ ^ ^ ^^^ 

Theorem 1: A network- wide scaling vector (3opt = (/?i^ j • ■ • ; /^w* ) ^^^^ maximizes the SNRt for a diamond network with 
the relay nodes performing ANC is given as 

'l3i,Tnax,i = k*, k* = a.rgma.x SNRt{l3it), 

Pi,max7 '^ t O , 

'j,max"'jt 



P\ 



opt 



1+ E Plmaxh'] 



^'^^ Jg^"^ ,• a Ok- 

^it 2-^ '^sj Pj^niax'^jt 

where S''* = {k*} U {i : /?f > A,™,„ i e AA \ {A:*}}- 

Based on our approach in this section to compute the optimal ANC rate in the Gaussian diamond networks, in the next 
section we introduce a greedy scheme to bound from below the maximum end-to-end rate achievable with analog network 
coding in the general layered networks. 

IV. General layered networks: a greedy scheme to lower bound the maximum ANC rate 

In a general layered network with L layers of relay nodes, consider layer 1,1 < I < L, and a node in the next I + V^ layer, 
denoted as f/+i or with a little abuse of notation as t. This scenario is depicted in Figure [3] For this subnetwork, for any 
scaling vector /3 we have 

SNRt = ^^(^^=i^'^'^»*)' (29) 

Using ([9| the problem of computing the maximum ANC rate for this subnetwork can be formulated as 

max SNRt, (30) 

where /3 = (/3i, . . . , /3n) and /3„a.T = {l^i^max ■■■, PN^max) with Pl^^^ = Pil¥.{siXs + ^^)^ i e N ,U = {I, . . . ,N}. 

Equating the first-order partial derivatives of the objective function with respect to Pi,i & JV, to zero, we get the following 
A^ + 1 conditions for local extrema: 

Y, S^P^hu =0 (31) 

/3^ = ^'^-^\^^> (32) 




► 



Fig. 3. A subnetwork of general layered network with L relay layers, depicting i* layer with N relay nodes and a node in the I + I'" layer. The received 
signal component at node i, 1 < i < Af, in the Z"' layer is denoted as SiX^, where x^ is the source symbol and corresponding noise component is denoted 



(signal component at destination t from all the nodes except node i) 



where 

jeM\{i} 

7i = 2_, Pj^^jt^i^i^j)/'^^^ (noise component at destination t from all the nodes except node i) 

jeAf\{i} 

As we established the extremal properties of conditions ( [T3| ) and ( [T4] l in Section III we can also prove that condition pT| ) 
leads to the global minimum of the objective function in ( pO] ) and the global maximum of the objective function occurs at 
one of the N hyperplanes defined by /3fc = Pk,max- Following a sequence of arguments similar to those used to establish 
Theorem n] for diamond networks, we can characterize the scaling vector for the nodes in the l^^ layer that optimally solve 
the problem ( pO] ) for the subnetwork under consideration. Note that in this subnetwork, the noises at different nodes in a relay 
layer are correlated, unlike the noises at relay nodes in the diamond network in Figure [2] This explains the difference between 
the SNR expression in ( p9l ) and the one in ( [TT| for the diamond network, and results in more complex analysis in the present 
case. 

Lemma 2: A scaling vector (3opt — {Pl^ ,■ ■ ■ , P'm ) ^^^ maximizes the SNRt for any subnetwork, as in Figure bl of the 
general layered network with the relay nodes in V^ layer performing analog network coding is given as 



/?; 



opt 






^i::^k\k* ^ argmax 5iVi?t(/3^), 

Sj + SiE{ J2 Zjl3j,maxhjt) - Ojjj 



hit{ajEzf/a^ 



^Ij) 



where fi'^ = (/3^', . . . ,;9^) with 



j,max ■} ' 



/3. 

s, + {siEzya 






^j^max'^jt 



a, 



7i 



huisjEzf/a-^ - s,E{z,Zj)/a^}pj 
2_\ SjPj^rnaxhjt, (signal component at destination t from the nodes in S^ ) 

/ I-' J, max"' J t 



f3j,maxhjt'E{ziZj)/a'^, (noise component at destination t from the nodes in S'^ ) 



and 5^-* = {k*} U {i : /3f > ft,„,,, i €M\{k*}}. 

Note that Lemma [2] reduces to Theorem [T] when the noise components at the relay nodes are uncorrected. 

Using Lemma 2 we can compute j3i '^ , the scaling vector for the nodes in the /* layer that maximizes the received SNR 
for node j, I < j < n;+i, in the I + I'*' layer. Among these n/+i scaling vectors for the nodes in the /* layer, let /Sj""" denote 



the one that solves the following problem 

jSl"'"' = argmax TT (1 + SNRk) (33) 



Pi, opt k=l 
l<j<ni 



The following corollary of Lemma 2 in LTJ establishes that among ni+i such scaling vectors, the scaling vector characterized 
by Pl""^ computes the tightest lower bound for the optimal value of the objective function in ( (33] l as well as ( [TO] i. 

Corollary 1 ([7], Lemma 2): Consider two scaling vectors /3; and 0i for the nodes in Z* layer If Hl-l^i (1 + SNRk)\„ > 
nfc=i'(l + SNRk)\^^, then SNRt{f3i) > SNRt{$i). 

Computing /S,'""" as above for each layer 1,1 < I < L,m conjunction with Corollary IT] allows us to construct a network-wide 
scaling vector /3/ou) = (/3i°™7 ■ ■ ■ , /3l" ) to compute a lower bouncrlto the optimal solution of (|9]). Formally, for a given layered 
network, f3iow is constructed as follows. 

Proposition 2: Consider a layered relay network of i + 2 layers, with the source s in layer '0', the destination t in layer 
'i + 1', and L layers of relay nodes between them. The Z* layer contains ni nodes, no = ^l+i = 1. A network- wide scaling 
vector f3iow — (/3i°'", • • ■ ,/3l"") that provides a lower bound to the optimal solution of (|9| for this network, can be computed 
recursively for 1 < / < L as 

(3'r - argmax [] (1 + SMRi+.^f^T , ■■■, /^t^Pltl')) 

/3;+j,-/ k=i 

l<j<ni 

Here /3j°"' = (/3j°™, . . . , (3l°^) is the vector of scaling factors for the nodes in the /* layer, and (3i 'f (computed using 
Lemma EJ) is the scaling vector for the nodes in the Z* layer that maximizes the received SNR for node j, I < j < ni+i, in 
layer I + 1. 

In the next section we analyze the performance of the greedy scheme of the Proposition [2] in the context of both a special 
class of layered networks and the general layered networks. 

V. Illustration 

We first demonstrate that the greedy scheme of Proposition l2] allows us to exactly compute the optimal ANC rate for a broad 
class of layered networks. Then, we give an example to show that for the general layered networks, the proposed scheme leads 
to the optimal rates within a constant gap from the cut-set upper bound asymptotically in the source power. 

Example 2 (A class of exactly solvable layered networks): Let us consider a class of symmetric layered networks where the 
channel gains along all outgoing links from a node are equal. An instance of such a network is obtained from the network 
in Figure HI when hgi = /is2,^i3 = hu, and /i23 — ^24. An implication of this property of the channel gains is that the 
received SNRs at every node in a layer are equal: SNRij ~ SNRi, I < j < ni,l < I < L. In this case, for each layer 
1,1 < I < L, fi\°^ computed in Proposition O^ is equal to the optimal (3°^ computed in [|7J Lemma 2]. Therefore, (3iow is the 
optimal solution of problem ( fTO] ) for this class of networks. 

Consider an instance of the network in Figure |4] when h^i = hs2, /113 ~ /ii4, and /i23 = Zi24- Such an instance belongs to 
the class of symmetric networks we are concerned with in this example. Using Proposition l2] the optimal solution of problem 
( [TOJ l for this instance is: 

h2Pl,maxhl 



Hopt 1 Pl.max, J c> 7_ , PZ,max, PA^vnax \, Wnere 



Pl.,raax ^2p^ _^ ^ 



S — ho{(3i^opthi + P2,opth2),Z = 1 + Pi^opthi + P2Mpt^2 

'clearly, choosing 0"°™ as in J33} for each layer / may lead, in general, to some performance loss at each layer as /3;°™ may not be the optimal vector 
of the scaling factors for the nodes in the layer I that solves 

"1+1 
argmax TT (1 + SNRk) 

The cumulative effect of this performance loss at each layer is that the end-to-end ANC rate computed at f3iow may not lead to the optimal solution of 
problem J9j. However, our results in the next section show that for a large class of layered networks there is no loss in the optimality and for other layered 
networks, the loss is small asymptotically in the network parameters. 




2 "24 ^ 

Fig. 4. General layered network with 2 relay layers between the source 's' and destination 't'. Each layer contains two relay nodes. 
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Fig. 5. Comparison of the ANC rate achievable with the scheme in Proposition l2] with the MAC upper bound for the layered network in Figure p] with 
Pi = P2 = Pg = P4 = 10, hi4 = /i24 = 2 and all other channel gains are equal to 10. Also plotted is the ANC rate when the scaling factors for all relay 
nodes are set to their respective upper-bounds. 



and we assume that Pihf > P2h\. ■ 

Example 3 (General layered networks): Let us consider the layered network of Figure |4] We compute a lower bound to 
the optimal ANC rate for this network using the greedy scheme in the Proposition l2] and compare it with the MAC upper 
bound in Figure [5] Also, plotted in this figure is the ANC rate achievable when the scaling factors for all relay nodes are set 
to their respective upper-bounds. We observe that in this case the ANC rate achieved with the greedy scheme of Proposition [2] 
approaches the capacity within one bit when Pg > 100. 

VI. Conclusion and Future Work 

We consider the problem of maximum rate achievable with analog network coding in general layered networks. Previously, 
this problem has been considered under certain assumptions on per node scaling factor and received SNR as without these 
assumptions the problem was presumed to be intractable. The key contribution of this work is a greedy scheme to exactly 
compute the optimal rates in a wider class of layered networks than those that can be addressed using prior approaches. In 
particular, using the proposed scheme for the Gaussian A^-relay diamond network, to the best of our knowledge, we provide the 
first exact characterization of the optimal rate achievable with analog network coding. Further, for general layered networks, 
our scheme allows us to compute optimal rates at most a constant gap away from the cut-set upper bound asymptotically in 
the source power. In the future, we plan to extend this work to non-layered networks, and to construct the optimal distributed 
relay schemes. 
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